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ABSTRACT 

We present optical interferometric polarimetry measurements of the Mira-like variables R Car 
and RR Sco, using the Sydney University Stellar Interferometer. By making visibility mea- 
surements in two perpendicular polarisations, the relatively low-surface brightness light scat- 
tered by atmospheric dust could be spatially separated from the bright Mira photospheric flux. 
This is the first reported successful use of long-baseline optical interferometric polarimetry. 
Observations were able to place constraints on the distribution of circumstellar material in 
R Car and RR Sco. The inner radius of dust formation for both stars was found to be less 
than 3 stellar radii: much closer than the expected innermost stable location for commonly- 
assumed astrophysical "dirty silicate" dust in these systems (silicate dust with a significant 
iron content). A model with the dust distributed over a shell which is geometrically thin com- 
pared to the stellar radius was preferred over an outflow. We propose dust components whose 
chemistry and opacity properties enable survival at these extreme inner radii. 
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1 INTRODUCTION 

The processes of dust formation and mass loss in evolved giants 
and supergiants remains poorly understood despite its manifold and 
profound implications to the cycle of matter through stars and back 
into the ISM. An extensive obs ervational and phenomenological 
literature (see e.g. lHabin3 for a review) bears witness to the 

barrage of astrophysical techniques applied to this problem, and 
progress has been considerable in many important areas. However, 
fundamental uncertainties remain over the mechanisms of eleva- 
tion of material, chemical nucleation of the dust, and subsequent 
outflow and global properties of the wind. As an example of the 
gap in our understanding, there has been to date no self-consistent 
physical model of the atmosphere of an evolved O-rich giant which 
encompasses the upper regions including even the most important 
basic phenomena: non-grey radiative transfer, pulsation, shocks, 
convection, non-equilibrium dust formation and the molecular at- 
mosphere. 

From an observational perspective, progress has been ham- 
pered until now by the difficulties inherent in imaging material on 
size scales of a few AU which encompass the photospheres them- 
selves and the dust formation and acceleration regions. Although 
modern long-baseline optical/IR interferometry is capable of ac- 
cessing these scales, imaging is not straightforward as targets can 
be complex, deeply embedded and often asymmetric. In particular 
for many Mira-type variables, the relatively low surface-brightness 
material in the outflow can be very difficult to detect against the 
high luminosity stellar flux. 

The stellar photosphere and circumstellar dust shell are usu- 



ally treated as completely separate regimes, as justified by the fol- 
lowing simple argument. As a first approximation, one might con- 
sider that dust grains in an optically thin envelope are grey or nearly 
so, and exist in thermal equilibrium with the radiation field. Fur- 
thermore, we can approximate the star as having a small angular 
extent, so that the mean intensity, J(A) at the dust condensation 
radius R c is equal to the luminosity of the central star L(X) di- 
vided by 16-K 2 Ilj:. In this case EquationQgives an estimate of the 
dust formation radius (here R s the radius of the central star, T a the 
effective temperature of the star and T c the dust condensation tem- 
perature). Given an effective temperature of 3000 K and a dust con- 
densation temperature of HOOK, the dust formation radius is 3.7 
stellar radii, which means that the physics of the photosphere and 
the circumstellar environment could be considered well-separated. 
Where dust is considered to be non-grey, but the dust formation ra- 
dius still lies well-outside the photosphere, conservation of energy 
and radiative equilibrium gives Equation|2| where Q is the absorp- 
tion coefficient of the dust and B is the Planck function in the form 
which gives power per unit area per steradian per unit wavelength. 
F urther discussion of these kinds of approximations can be found 
in lSchutte & Tielenj Jl989l) . 



= / /Q abs (A)L(A)dA 
c y 1671-2/ Q abs (A)B(A,T c )dA 

The absorption coefficients for dust around O-rich Mira are 
often assumed to be those of "dirty" silicates I lones & Merrill 
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ll976tlOssenkopf et alll993) . These silicates have optical constants 
very similar to am orphous olivine Mg^Fei-^SiC^ with x — 0.5 
JPorschner et alJll995h . with the exception of the spectral range 
2-8 /im, where inclusions such as solid Fe are required to match 
observationally derived optical constants. Olivine with x = 0.5 
can form only slightly closer to the star than grey dust. Using the 
emitted flux from a 3000 K effective temperature Mira model from 
Ireland et al. (2004a) (the M09n model), the optical constants from 
iDorschner et alj 1 19951) and an assumed condensation temperature 
of 1 100 K, this olivine can form at a radius of 6.4 stellar radii: still 
a comfortable distance from the upper atmosphere of the star. 

Apart from the optical properties, the other factor governing 
grain survival at smaller radii is the sublimation temperature. The 
dust species often considered to be the first to form in higher effec- 
tive temperature stars from semi-empirical considerations is corun- 
dum lEga n&Sloarl200ll) . AI2O3 , which is stab le at considerably 
higher temperatures than silicates lSalpeterll977l) . Its formation ra- 
dius for this same 3000 K model s tar is 2 . 1 stellar radii, using the 
optical constants from Koi ke et all 1 19951) and an assumed conden- 
sation temperature of 1400K. 

The multi-wavelength diameter measurements of six Mira- 
like variables made by Ireland et al. (2004b) show an increase in 
apparent diameter from 900 nm to 700 nm with a component which 
appears uncorrelated with the strength of TiO bands. These obser- 
vations are best explained by scattering from dust closer than a few 
stellar radii (although molecular line blanketing makes quantitative 
interpretation difficult). The presence of extensive dust this close 
is, as discussed above, difficult to reconcile with the conventional 
wisdom on dust formation around AGB stars. 

Despite the fact that the net polarisation arising from a spher- 
ically symmetric dust shell is zero, it is possible to separate it from 
the unpolarised photosphere with a polarisation-sensitive high- 
resolution measurement. Scattering from the dust shell will have 
a polarisation at any point with the E field in an azimuthal direction 
with respect to the stellar centre. This is due to the simple relation- 
ship that scattered light from small particles in the Rayleigh limit 
is polarised orthogonal to the plane containing the source, scatterer 
and observer. 

If we now consider the response of a long-baseline interfer- 
ometer which measures the Fourier components (called visibilities 
when total flux is normalised to unity) of such an image, these will 
have different amplitudes and phases in different polarisations. For 
a partially-resolved circularly symmetric source, the linear polari- 
sation with the E field parallel to the baseline of observation will 
have a lower amplitude and the same phase as the visibility mea- 
sured using the polarisation perpendicular to the baseline. An at- 
tempt was made using this technique to detect scattering by free 
electrons around j3 Orionis using the long-baseline Narrabri Stel- 
lar Intensity Interferometer iBrown et alll974l) . but the instrument 
lacked the required sensitivity. Here we report the first successful 
use of this technique on a long-baseline optical interferometer. The 
observational technique is described in the following section, while 
Section[3]gives the results and Section|4]a discussion and astrophys- 
ical interpretation. 



2 OBSERVATIONS 

The Sydney University Stellar Interferometer (SUSI) is a long 
baseline optical interferometer with presently operational baselines 
from 5 to 160 m (further details can be found inD avis et alJ ll9991). 
The observations presented here used 5 and 10 m baselines and 



a filter with central wavelength of 900 nm and full-width half- 
maximum 80 nm. This filter was chosen so as to minimise con- 
tamination from the strongest features caused by the TiO molecule. 
It is in principle possible to measure visibilities at arbitrary polari- 
sations through the use of quarter-wave plates and polarisers at ar- 
bitrary orientations. For a discussion of a general way to calibrate 
a long b aseline interf erometer for imaging in all Stokes parame- 
ters, see Elias (2004). Due to polarisation-dependent phase-shifts 
and mirror reflectivities over the many mirrors in the SUSI opti- 
cal chain, arbitrary initial polarisation states are not well preserved. 
The measurements therefore are restricted to states with pure ver- 
tical and horizontal electric field orientation with respect to the op- 
tical path. While they will be attenuated, these states should not 
suffer mixing from any reflection with the exception of the initial 
reflection off the siderostat when observing away from the merid- 
ian. Care was taken to only observe stellar targets when near the 
meridian. Given SUSI's North-South baseline orientation, this ob- 
serving strategy resulted in visibility measurements with the E field 
at all times being within 19 degrees of perpendicular (Vx) or par- 
allel (Vj|) to the baseline. For crucial observations where a direct 
comparison was made between the visibilities in these two polari- 
sations, data were restricted further to lie within 8 degrees of these 
orientations. 

The beam combining system used for these observations con- 
sisted of a pupil-plane beam combiner feeding two Avalanche 
Photo-Diode (APD) detectors. By repeatedly scanning the optical 
delay linearly through the white-light fringe position, the light de- 
tected from each output of the combining beamsplitter is modulated 
at a frequency that depends on the scanning rate and the observ- 
ing wavelength, with the fringe signals 180 degrees out of phase 
at the two output s . A d escription of a similar system can be found 
in Bald win et al] ll994h . and further details of the beam combin- 
ing sy stem used for these observations will be found in lDavis et alJ 
(2005) (in preparation). Although the two outputs are usually used 
differentially to reduce the effects of scintillation noise, it is pos- 
sible to obtain an estimate for the fringe visibility from only one 
output, as in FigureQ The use of only one detector means that the 
bias in the power spectrum is not flat due to scintillation, and in- 
correct subtraction of this bias becomes an additional error source 
for low V 2 . The observational strategy employed consisted of two 
separate parts: the first aimed at measuring V± and Vji separately, 
while the second measured the ratio Vx/Vu. These are described in 
turn below. 

Firstly, polarisers were co-aligned and placed in front of the 
two outputs of the beam combiner, enabling a high signal-to-noise 
measurement of the correlation in linear polarised light. As is usual 
with optical interferometry measurements, observations of targets 
were interleaved with un- or partially-resolved reference stars to 
calibrate the system correlation (the V 2 response to a point source) 
as a function of time. Before use in adjusting the science mea- 
surements, this system response was also corrected for the known 
nonzero diameters of the reference stars given in Table Q In this 
phase of the observational program, measurements were predomi- 
nantly made of Vx which had a much higher signal-to-noise than 
Vji . This is mainly due to partially polarising beamsplitters used at 
SUSI to split off light to the tip/tilt adaptive optics camera, which 
had the effect of significantly lowering the system throughput for 

The second part of the observing strategy consisted of simul- 
taneous observations of Vx and Vji , made possible by placing po- 
larisers with orthogonal orientations in front of the two outputs of 
the beam combiner. These observations gave a measurement of the 
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Table 1. Assumed diameters for the calibrator stars used for these observa- 
tions 



Star 


Uniform Disk (UD) 


Method 2 




Diameter (mas) 




HR6241 


5.9 


VB99 


HR 6630 


4.2 


VB99 


HR 4050 


4.9 


C99 


HR 6553 


2.2 1 


082 


HR 3685 


1.54 


082 


HR 3884 


2.5,2.8 


S05 



1 This star also has a companion 3.39 magnitudes fainter and at distance 
of 6.47" (figures taken from the H1PPARCOS catalogue). This was over- 
resolved by SUSI and taken into account in the calibrated values f or V 2 . 

2 VB99: Diameters estimated from the B-K relat ionship fromlvan Belld 
i 1999), using the photometry from the catalogue of Ducati 1 2002). 

C99: Diameter taken from Cohen et al. 1 1999), with limb-darkening correc- 
tion from lClareJ feOOCl) , 

082: Diameters taken from lOchsenbein & HalbwachJ Il982t) . 
S05: Diameters of this Cepheid taken from near-simultaneous measure- 
ments at SUSI, at early April and early May epochs, consistent with the 
published diameters of Kervella et al. l2004> . 

ratio V±/V\\, and due to the simultaneous data, the ratio could be 
obtained independent of seeing effects so that rapidly interspersing 
calibrator star observations was no longer necessary, and longer 
total integration times were possible. Calibrator stars, observed at 
a different time of night, were still used for these observations to 
correct for instrumental effects that might affect this visibility ra- 
tio. Some calibrators were moderately resolved, and calibrators of 
a range of spectral types were used, however there was no evidence 
for systematic variation of the visibility ratio Vj_/Vji between cali- 
brator stars. 

A typical observation in a single polarisation consisted of 
4000 50 ms sweeps through the fringe envelope, while the simul- 
taneous observations of V± and Vji contained up to 20000 50 ms 
sweeps through the fringe envelope. The slew and acquisition time 
between source and calibrator generally averaged 3-4 minutes. 
Details of all observations can be found in Table [2] Note that 
squared visibility (or correlation) quantities are listed, because in 
the photon-noise limited regime, squared visibility must be aver- 
aged over many scans and error distributions are closest to Gaus- 
sian in these quantities. 

There are several sources of systematic error in the measure- 
ment of V { f/Vl. One error source which was corrected for in the 
values in Table [2] was a 5nm shift in effective centre wavelength 
between the filter in each output channel. The effect of this was 
to make the souce appear slightly more resolved in one channel 
then the other. This effect was well characterised by observing in 
a mode where both output channels observed the same polarisation 
state, and resulted in a maximum correction of 3% in Vjj 2 / Vj 2 for 
the most well-resolved observation. 

Other errors, not included as their effects were computed to be 
small, consisted of imperfect polariser alignment, offsets between 
the actual polarisation passed and the ideal desired sky orientation, 
and depolarisation and mixing due to alignment errors in the beam 
train. With the assumptions of spherical symmetry and rayleigh 
scattering, each of these effects multiplies Vji — VI by a factor 
7 which is has a modulus smaller than 1. We will define <f> to be 
the azimuthal coordinate describing the rotation of the siderostat 
normal about a horizontal north-south axis, with <f> = describing 
the meridian. The effect of vertical/horizontal linear polarisation 
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Figure 1. An example power spectrum (for the calibrator star HR 6241 in 
the 900 nm filter) obtained from averaging the power spectra from 3000 
fringe scans together. Power is in arbitrary units. Note that although the 
scintillation power is clear at low frequencies, bias subtraction is relatively 
simple, and in this case a constant bias (dashed line) would be a good ap- 
proximation over the range of frequencies covered by the fringe peak. 

rotation with respect to the baseline orientation then gives a fac- 
tor 71 = cos(2(j!>). The mixing of linear with circular polarisation 
states due to phase changes on reflection gives a second factor 72 . 
If we define S to be the phase shift between s and p polarisations at 
the siderostat mirror (itself a function of the angle of incidence of 
the starlight), then 72 is given by: 

72 = cos (5/2) + sin 2 (<5/2) cos(4</>) (3) 

For example, S is 37 degrees for SUSFs siderostats at 900 nm 
and an angle of incidence of 45 degrees, which gives 7 = 0.991 for 
the maximum of 6 degrees for dual-polarisation measurements 
RR Sco. Taking all of these effects into account, we estimate that 
for no observation will the true values for V^/Vj 3 be more than 
0.6% higher than the measured values given in Table H This is 
smaller than statistical measurement errors. 



3 RESULTS AND MODEL FITTING 

As a first step in interpretation of the results, from the observed 
quantities in Table|2|we derive the quantities Vji — V± and V±. For 
R Car, we have averaged data into four points at projected base- 
lines of 4.2 and 8.4 m and covering pulsational phases of ~ 0.08 
and ~ 0.15. Similarly, observations of RR Sco have been distilled 
into measurements at baselines of 5.0 and 10.0 m with a phase of 
0.95. As the signal-to-noise of the initial observations is high, the 
transformation from the observed quantities in V 2 to Vj| — V± and 
V± could be performed while preserving the symmetric Gaussian 
error distributions reasonably well. 

We assume that thermal photons from the star can be well rep- 
resented by a Uniform Disk (UD) visibility curve, and that the dust 
distribution is spherically symmetric. The first assumption that the 
star is UD-like is reasonable because at phases near maximum, the 
atmosphere of a Mira variable is expected to be relatively com- 
pact (e.g. see Ireland et al. (2004a)), and the filter was chosen so 
that there was minimal contamination by the TiO molecule. The 
second assumption that the dust is spherically symmetric is reason- 
able for R Car because of its low total polarisation. Although any 
measured polarisation may be interstellar, a low total polarisation 
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Figure 2. Fits to the first epoch (phase 0.08) of R Car (upper plots), the second epoch (phase 0.15) of R Car (middle plots) and the single epoch (phase 0.95) 
of RR Sco (lower plots) with data reduced to Vj| — V± (left plots) and V± (right plots). The solid line is the thin shell model fit, the dashed line is the outflow 
model fit and the dotted line is the uniform disk only fit. 
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Table 2. Summary of SUSI observations of R Car and RR Sco. Visual phase (Col 4) were obtained from AAVSO data. The projected baseline B p and position 
angle range are given in Cols 5 & 6, while the type of measurement and its recorded value are in Cols 7 & 8 



Date Star JD Phase B p (m) PA range (°) Measurement Value 



2004 Apr 8 


RCar 


2453104 


0.08 


4.2 


-15. .-6 


vl 


0.582 ± 0.04 


2004 Apr 8 


RCar 


2453104 


0.08 


8.4 


1..18 


vl 


0.123 ±0.015 


2004 Apr 8 


RCar 


2453104 


0.08 


4.2 


-5..0 


V,? A/2 


1.160 ±0.03 


2004 Apr 9 


RCar 


2453105 


0.08 


4.2 


-10..-2 


vl 


0.556 ±0.04 


2004 Apr 9 


RCar 


2453105 


0.08 


8.4 


5.. 10 


vl 


0.125 ±0.015 


2004 Apr 9 


RCar 


2453105 


0.08 


8.4 


-1..4 


v l\/ v l 


1.220 ±0.06 


2004 May 1 


RCar 


2453127 


0.15 


4.2 


-12..0 


vl 


0.592 ±0.04 


2004 May 1 


RCar 


2453127 


0.15 


4.2 


3. .8 


v?/vl 


1.095 ±0.03 


2004 May 2 


RCar 


2453128 


0.15 


8.4 


-4..0 


v l 


0.129 ±0.01 


2004 May 2 


RCar 


2453128 


0.15 


8.4 


-4..-1 


vpvl 


1.150 ±0.11 


2004 May 6 


RCar 


2453132 


0.16 


8.4 


0..5 


v \\/ v l 


1.220 ±0.10 


2004 Jul 28 


RRSco 


2453215 


0.95 


10.0 


4..7 


v t 


0.196 ±0.012 


2004 Jul 28 


RRSco 


2453215 


0.95 


10.0 


-10..-7 




0.259 ±0.015 


2004 Jul 28 


RRSco 


2453215 


0.95 


10.0 


-5. .2 


v^/vl 


1.195 ±0.05 


2004 Jul 29 


RRSco 


2453215 


0.95 


5.0 


-11. .10 


vl 


0.646 ± 0.01 


2004 Jul 29 


RRSco 


2453215 


0.95 


5.0 


-2..5 


v \\/ v l 


1.063 ±0.014 



signal places a rough upper-limit on the amount of light scattered 
by asymmetrical regions close st to the star. Measured polarisation 
fractions are 0.62% at 652 nm jMcLean & Clarke! 197/1) . 0.52% at 
652.5 nm ICodina-Landaberrv & Maealhaesl 19801) and 0.5 1 % in V 
band ISerkowski & Shawl20oHT with all authors measuring a sig- 
nificant downward trend in polarised fraction with increasing wave- 
length. These numbers can be compared with the roughly 15% total 
scattered light as measured in this paper at 900 nm. No such com- 
prehensive polarimetry measurements could be found in the litera- 
ture for RR Sco. 

Two dust models were fit to the data: a thin shell and an out- 
flow. For the thin shell model dust is assumed to exist in a static 
infinitesimally thin shell, while for the outflow model dust forms at 
some condensation radius and is then radiatively accelerated pro- 
ducing a stellar wind. For the outflow model, the dust density as a 
function of radius can be derived as (Schutte & Tielens 1989): 

p(r)otl/{(r-(l-v 2 )) 1/2 r S/2 } (4) 

Here r is radius in units of the condensation radius and v c is 
the outflow velocity at r c as a fraction of the final outflow velocity. 
We only consider the limiting case of v c = in this paper as it 
is provides a better fit to the data in all cases. After fixing v c to 
zero, both models have only two free parameters, the radius of dust 
formation and the optical depth of the dust shell. In both models, we 
assume that the dust grain radii are much smaller than the observing 
wavelength, so that the Rayleigh scattering approximation can be 
used, and the light scattered by the dust is all assumed to originate 
from the central star (i.e. only a single scattering is included). 

The polarisation signal from a dust shell illuminated by an 
extended photosphere will be weaker than that illuminated by a 
compact object (or point source) of equivalent luminosity at the 
stellar centre. We proceed here to derive the dilution of the signal 
due to the geometrical extension of the underlying star. Consider a 
dust grain at a distance from the centre of a star whose radius 
is r 3 . The polarisation properties of light scattered from this grain 
will vary according to the angle of the incident photon, which in 
turn depends on the original location of the emission on the surface 
of the star. The net polarisation is then found by integration over 
the disk of the star visible to the dust grain. 




Figure 3. An example of the geometry for which the polarised intensity was 
calculated, using the approximation of Rayleigh scattering from an optically 
thin shell (see text). 

For Rayleigh scattering, if the angle between the polarisation 
vector of the radiation seen by the observer and the ray incident on 
the dust grain is a, then the intensity of the scattered radiation is 
proportional to sin 2 a. Let 8 be the angle between the ray incident 
on the dust grain and the direction from the dust grain to the star's 
centre, let cf> be an azimuthal coordinate describing a ring on the 
stellar surface and let f3 be the angle between the vector from the 
star's centre to the dust grain and the vector from the dust grain to 
the observer. An example of this geometry when all vectors are co- 
planar is shown in Figure|3| Then, using spherical trigonometry, a 
can be expressed in terms of 9, <f), and /?. The integral that gives the 
scattered intensity is: 

/•2-7T />arcsin — ^> 

1=1 d(f> I sin 6 sin 2 add (5) 

Jo Jo 
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We consider the polarised intensities of scattered light from 
the dust grain decomposed into those with electric fields in radial 
(I r ) and azimuthal (I a ) directions with respect to the centre of the 
underlying disk. Evaluating the integral for I r and I a gives: 



It 
Ia 

where u 



2 — u — u 2 , u + u 2 , 



■ cos P 



4 + u + u 2 
6 



(6) 
(7) 
(8) 



These integrals are normalised so that forward scattering from 
a point source gives intensities of 1. Note that for dust distant from 
the star, u — 1 and this formula reduces to the familiar formula 
for Rayleigh scattering. These intensities were numerically inte- 
grated for the two distributions of dust discussed above to give a 
two-dimensional image, and a two-dimensional Fourier transform 
applied to this image to give the visibility functions arising from 
the dust scattering. Radiation back-scattered onto the surface of the 
star was assumed to be absorbed. After normalising these dust visi- 
bility functions to 1 at zero baseline, the overall visibility functions 
become: 



Vii 



V, 



V s + {l-x b )(l 



1 - x b (l - e-T*) 
e- T «V s + {I - x b )(l - e 



<)Vd± 



1 -351,(1 • 



') 



(9) 



(10) 



Here Td is the optical-depth of the dust, Vd±, Vm and V s 
are the visibility functions for the dust in both polarisations and 
the star respectively, and Xb is the fraction of the scattered light 
backscattered to re-intercept the surface of the star. The derived 
fit parameters for both stars are given in Table [5] with the fits for 
both stars shown graphically in Figure [5] In addition to quantities 
derived from fitting to the high-resolution data, for comparison in 
Table|5|we also present expected stellar diameters based on the ob- 
served luminosities. These are given in the Model D 3 column, and 
were derived using near-maximum dust-free M series atmospheric 
models from Ireland et al. (2004a) (interpolating between the M08, 
M09n and M10 models), with the distance of the model star set to 
match the K band magnitudes measured by Whit elock et alJ fcOOtl) 
at the same phas es. Note that the J-K colours of these stars as mea- 
sured by Whitelo ck et alJ 12000) are consistent with this model se- 
ries within the observational cycle-to-cycle scatter. 



4 DISCUSSION 

The dust free (UD only) fits in Table [5] are clearly eliminated for 
both stars. The high % 2 values have their origin in the inability of 
this scattering-free model to generate any polarisation signal. Of 
the two models with dust, the thin shell model is a better fit than 
the outflow model for both stars. The best-fit condensation radii in 
the outflow models are also uncomfortably close to the continuum- 
forming photospheres for both stars. Although there may be dust 
types that could plausibly exist in thermal equilibrium with the radi- 
ation field at these radii (as discussed below), the gas temperatures 
in the model star of Section 1 are 2000 K at 1.1 stellar radii , much 
too hig h for grain nucleation (e.g. see the discussion in ljeong et alJ 
(2003)). The only way dust could exist so close to the continuum- 
forming photosphere is if it had formed at near-minimum phases 



and then fallen in. However, this kind of motion is in conflict with 
the general assumptions of the outflow model. 

The two epochs of R Car do not show statistically significant 
differences, but the small differences at the level of about l-cr are 
in the expected direction. The star appears to increase in apparent 
size from phase 0.08 to 0.15, as predicted by the models, and the 
optical-depth in scattering appears to decrease, consistent with the 
dust being part of an infalling layer that is partially sublimating at 
these near-maximum phases. For RR Sco, the discrepancy between 
the luminosity/model-predicted diameter and the best fit diameter 
is almost certainly due to contamination from the TiO molecule. 
A moderate strength TiO band at 850 nm is included in the wing 
of the SUSI filter profile which could account for the enlargement 
observed, although predicting the effects of such molecular con- 
tamination is made difficult by its strongly cycle-dependent nature 
(Ireland et al. 2004a). 

From the discussion of dust formation radii in Section 1, it 
is clear that "dirty" silicates cannot form a major fraction of the 
optical depth in scattering due to dust around R Car. For survival 
within the hostile environment close to the photosphere, dust must 
have very low absorption between 1 and 4 fim where most of the ra- 
diation from the central star is emitted, and higher absorption coef- 
ficients at longer wavelengths promoting efficient radiative cooling. 
Corundum is a clear possibility here, but another obvious dust type 
that satisfies these criteria is forsterite, the Mg-rich olivine that is 
predicted to be the first significant silicate to form in chemical equi- 
librium calculations l Gail & Sedlmavr 1999). The optical constants 
for this dust species as given in |ja^e^talTHo03l) result in absorp- 
tion coefficients 300 times higher at 10 /im than at 2 /im. This dust 
is stable at a radius of 1.2 stellar radii from the 3000 K model star 
of Section 1. As a plausible and abundant dust candidate, we will 
assume that the dust around these stars is made of forsterite for the 
following analysis of the implications of the dust model parame- 
ters. 

Low absorption at wavelengths where the bulk of the stellar 
radiation is emitted also results in lower radiation pressure on the 
dust. Assuming full dynamical coupling between the gas and dust 
and full Mg condensation, we can calculate the radiative accelera- 
tion as a fraction of the gravitational acceleration for an optically- 
thin model dust shell. For a shell made of forsterite around the 
3000 K model star of Section 1 using Mie theory to calculate ab- 
sorption and scattering coefficients, this fraction is 41% for spheri- 
cal grains of radius 0.1 /im and 1.5% for grains of radius 0.02 /im. 
We therefore propose that this dust is formed in material primarily 
elevated by shocks, which in turn means that significant temporal 
variation in the dust shell optical depth and/or radius is expected. 

Where the Rayleigh limit of small particles with respect to the 
wavelength of scattered light applies, it is sufficient to approximate 
dust with a distribution of sizes as a homogeneous population with 
a single effective radius a e . For dust composed of forsterite, with 
a scattering optical depth of 0.15 at 900 nm, typical of the obser- 
vations presented here, the optical depth in absorption at 10 /im 
is 1 for a e = 55 nm. Results from a radiative transfer model of a 
geometrically thin shell of forsterite with a fixed temperature of 
1 100 K and an optical-depth of 1 at two stellar radii from the centre 
of a spherical black-body 3000 K star are displayed in Figure|4| It 
can be seen that there is no strong mid-infrared feature evident in 
either emission or absorption due to the dust shell. It would seem 
that the dust here represents an intermediate case between the com- 
pact geometries typically yielding absorption features, and spectra 
exhibiting emission which is usually associated with extended ge- 
ometries. This result was found to be robust, and insensitive to the 



Dust Scattering in Miras R Car and RR Sco 7 



Table 3. Fit parameters for both epochs of R Car and the single epoch of RR Sco, giving the apparent stellar diameter D s , a model value for D s (see text), the 
diameter of the dust condensation shell D c the optical depth at 900 nm T900 and the reduced \ 2 f° r me fit- Note that there is only one degree of freedom for 
the Thin Shell and Outflow models. Details of models given in the text. 



Star 


Phase 


Fit Type 


D s 


Model D s 


D c 


T"900 


x 2 








(mas) 


(mas) 


(mas) 






RCar 


0.08 


Thin Shell 


15.8 ±0.6 


14.9 


32.3 ± 1.9 


0.19 ±0.03 


1.2 


RCar 


0.08 


Outflow 


15.6 ±0.6 


14.9 


17.3 ±3.1 


0.38 ± 0.07 


4.5 


RCar 


0.08 


UD Only 


18.7 ±0.3 


14.9 






15.4 


RCar 


0.15 


Thin Shell 


16.6 ±0.6 


16.0 


31.3 ±3.6 


0.14 ±0.04 


0.1 


RCar 


0.15 


Outflow 


16.2 ± 0.6 


16.0 


18.3 ±3.8 


0.26 ± 0.06 


0.9 


RCar 


0.15 


UD Only 


18.4 ±0.3 


16.0 






6.3 


RRSco 


0.95 


Thin Shell 


12.3 ±0.4 


8.3 


24.0 ±2.0 


0.14 ±0.03 


1.1 


RRSco 


0.95 


Outflow 


12.3 ±0.4 


8.3 


13.5 ± 2.2 


0.22 ± 0.04 


3.7 


RRSco 


0.95 


UD Only 


14.5 ± 0.2 


8.3 
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optical-depth of the shell. Note, however, that a shell of forsterite 
with a 10 micron optical-depth much greater than 1 would heat up 
beyond its sublimation temperature at the inner edge of the thin 
dust shell. In order to meet this criteria while maintaining the ob- 
served 900 nm optical depth, the effective particle size a e is con- 
strained to be larger than about 50 nm. From the point of view of 
energetics, the failure of the dust to have a profound impact on 
the spectrum is simply a consequence of its optical properties: its 
inability to absorb in the near-infrared means that it cannot redis- 
tribute the bulk of the stellar flux. The dust scattering measured by 
the technique in this paper is thus largely independent of the form 
of the infrared spectrum, and could therefore have remained hid- 
den from investigations based on SED fitting. By this argument, 
the mid-infrared excess in the spectrum must come from a mate- 
rial that absorbs more strongly than forsterite in the near-infrared. 
Dust that has evolved and been enriched with iron (that condenses 
at lower temperatures) in an outflow consisting of shells partially 
ejected in previous pulsation cycles could satisfy this property. 
More detailed mid-infrared spectral modeling is beyond the scope 
of this paper, but would need to include the emission from the ex- 
tended molecular atmosphere as predicted by self-consistent mod- 
eling by I Jacob & Scholzl 12 002) and semi-empirical modeling of 
IWemerli2004 and lOhnakj 120041) . 

The dust uncovered here is not incons istent with the mid- 
infrared measurements of Dan chi et alj Il994t) . These authors fit 
an outflow model to R Leo (which has slightly more pronounced 
mid-infrared emission than either R Car or RR Sco), with an op- 
tical depth at llpim of 0.1. For forsterite dust with a c = 55 nm, 
this corresponds to a 900 nm optical depth of 0.025. If a similar 
outflow were to exist around RR Sco or R Car in addition to the 
shell model, then its effect on the fit parameters in Table [5] would 
be roughly within the quoted errors. In order to measure scattering 
from such an outflow, shorter baselines and higher precision would 
be required. 

Both R Car and RR Sco have rel atively low mass-loss rates 
as constrained by CO line observations. Groenew egen et all l 1999) 
modeled the mass loss rate for RR Sco to be 1.1 x 10 -8 Mq per 
year and that of R Car to be less than 1.6 x 10~ 9 Mq per year. 
If we assume full Mg condensation at solar metallicity, then our 
best fit model with a shell at 2 stellar radii and an optical-depth 
0.15 in scattering at 900 nm has a total mass of 2.8 x 10~ 6 M Q if 
we continue our assumptions of dust composed of forsterite with 
a e = 0.55 nm and full condensation. Less than full condensation, 
or dust composed of the lower-abundance corundum would give 
an even higher mass shell. This means that the observed mass-loss 
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10 15 20 25 
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Figure 4. An example mid-IR spectrum (solid line) for a 1 100 K thin-shell 
model with a 10 fim optical-depth of 1 at 2 stellar radii, illuminated by a 
star with a 3000 K black-body spectrum. Flux is normalised to the 8 fim 
flux, and the 3000 K black-body spectrum is overlaid for reference (dotted 
line). Normalised IRAS Low Resolution Spectrometer (LRS) spectra are 
overplotted for R Car (crosses) and RR Sco (diamonds). The emission fea- 
tures in these spectra almost certainly come from a region geometrically 
more distant than that probed here (see text). 



rates are only consistent with the thin-shell model as long as the 
shell is considered nearly static, taking hundreds or thousands of 
pulsational cycles for material in the shell to be fully ejected. This 
is in turn consistent with the low radiative acceleration expected 
from dust that can exist so close to the stars. 

Outflow models, on the other hand, were not so easy to fit 
within a self-consistent picture of stellar mass loss. Using the mea- 
sured CO outflow velocity f rom l Yourigl 1 1 995|) of 3 km/s for RR Sco 
we can calculate the mass-loss rate from a typical outflow model in 
Table|3| Again assuming full Mg condensation at solar metallicity 
for forsterite dust with a c = 55 nm, the mass loss rate for a model 
with our best fit optical-depth 0.25 in scattering at 900 nm and con- 
densation radius of 1.1 stellar radii is 2.7 x 10~ 7 Mq per year, 
much higher than the observed mass-loss rates. This demonstrates 
once again that the outflow model can be ruled out. 
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5 CONCLUSION 

Using optical interferometric polarimetry, we have spatially sepa- 
rated the component of flux at 900 nm scattered by dust around the 
Mira variables R Car and RR Sco from their photospheric emission. 
We found that the inner radius of dust formation around these stars 
to be less than three stellar radii, consistent with dust that is rela- 
tively transparent between 1 and 4 fim, such as iron-poor silicates 
or corundum. This dust exists in a shell-like structure around these 
stars that may have little influence on their mid-infrared spectra and 
is not part of an outflow. This result demonstrates the complexity 
of the circumstellar environment of Mira variables where mass loss 
and dust formation are only made possible by pulsation. Simple 
outflow models are grossly inadequate, and it is necessary to con- 
sider the changing optical properties of dust as it evolves from the 
condensation radius outwards. 
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